A national-scale assessment of land subsidence in China's major cities.

China's massive wave of urbanization may be threatened by land subsidence. Using a spaceborne synthetic aperture radar interferometry technique, we provided a systematic assessment of land subsidence in all of China's major cities from 2015 to 2022. Of the examined urban lands, 45% are subsiding faster than 3 millimeters per year, and 16% are subsiding faster than 10 millimeters per year, affecting 29 and 7% of the urban population, respectively. The subsidence appears to be associated with a range of factors such as groundwater withdrawal and the weight of buildings. By 2120, 22 to 26% of China's coastal lands will have a relative elevation lower than sea level, hosting 9 to 11% of the coastal population, because of the combined effect of city subsidence and sea-level rise. Our results underscore the necessity of enhancing protective measures to mitigate potential damages from subsidence.

Carbon fluxes of China's coastal wetlands and impacts of reclamation and restoration.

Coastal wetlands play an important role in regulating atmospheric carbon dioxide (CO2) concentrations and contribute significantly to climate change mitigation. However, climate change, reclamation, and restoration have been causing substantial changes in coastal wetland areas and carbon exchange in China during recent decades. Here we compiled a carbon flux database consisting of 15 coastal wetland sites to assess the magnitude, patterns, and drivers of carbon fluxes and to compare fluxes among contrasting natural, disturbed, and restored wetlands. The natural coastal wetlands have the average net ecosystem exchange of CO2 (NEE) of -577 g C m-2 year-1, with -821 g C m-2 year-1 for mangrove forests and -430 g C m-2 year-1 for salt marshes. There are pronounced latitudinal patterns for carbon dioxide exchange of natural coastal wetlands: NEE increased whereas gross primary production (GPP) and respiration of ecosystem decreased with increasing latitude. Distinct environmental factors drive annual variations of GPP between mangroves and salt marshes; temperature was the dominant controlling factor in salt marshes, while temperature, precipitation, and solar radiation were co-dominant in mangroves. Meanwhile, both anthropogenic reclamation and restoration had substantial effects on coastal wetland carbon fluxes, and the effect of the anthropogenic perturbation in mangroves was more extensive than that in salt marshes. Furthermore, from 1980 to 2020, anthropogenic reclamation of China's coastal wetlands caused a carbon loss of ~3720 Gg C, while the mangrove restoration project during the period of 2021-2025 may switch restored coastal wetlands from a carbon source to carbon sink with a net carbon gain of 73 Gg C. The comparison of carbon fluxes among these coastal wetlands can improve our understanding of how anthropogenic perturbation can affect the potentials of coastal blue carbon in China, which has implications for informing conservation and restoration strategies and efforts of coastal wetlands.

Air quality improvements can strengthen China's food security.

Air pollution exerts crucial influence on crop yields and impacts regional and global food supplies. Here we employ a statistical model using satellite-based observations and flexible functional forms to analyse the synergistic effects of reductions in ozone and aerosols on China's food security. The model consistently shows that ozone is detrimental to crops, whereas aerosol has variable effects. China's maize, rice and wheat yields are projected to increase by 7.84%, 4.10% and 3.43%, respectively, upon reaching two air quality targets (60 µg m-3 for peak-season ozone and 35 µg m-3 for annual fine particulate matter). Average calories produced from these crops would surge by 4.51%, potentially allowing China to attain grain self-sufficiency 2 years earlier than previously estimated. These results show that ozone pollution control should be a high priority to increase staple crop edible calories, and future stringent air pollution regulations would enhance China's food security.

Increased photosynthesis during spring drought in energy-limited ecosystems.

Drought is often thought to reduce ecosystem photosynthesis. However, theory suggests there is potential for increased photosynthesis during meteorological drought, especially in energy-limited ecosystems. Here, we examine the response of photosynthesis (gross primary productivity, GPP) to meteorological drought across the water-energy limitation spectrum. We find a consistent increase in eddy covariance GPP during spring drought in energy-limited ecosystems (83% of the energy-limited sites). Half of spring GPP sensitivity to precipitation was predicted solely from the wetness index (R2 = 0.47, p < 0.001), with weaker relationships in summer and fall. Our results suggest GPP increases during spring drought for 55% of vegetated Northern Hemisphere lands ( >30° N). We then compare these results to terrestrial biosphere model outputs and remote sensing products. In contrast to trends detected in eddy covariance data, model mean GPP always declined under spring precipitation deficits after controlling for air temperature and light availability. While remote sensing products captured the observed negative spring GPP sensitivity in energy-limited ecosystems, terrestrial biosphere models proved insufficiently sensitive to spring precipitation deficits.

Non-symmetric responses of leaf onset date to natural warming and cooling in northern ecosystems.

The northern hemisphere has experienced regional cooling, especially during the global warming hiatus (1998-2012) due to ocean energy redistribution. However, the lack of studies about the natural cooling effects hampers our understanding of vegetation responses to climate change. Using 15,125 ground phenological time series at 3,620 sites since the 1950s and 31-year satellite greenness observations (1982-2012) covering the warming hiatus period, we show a stronger response of leaf onset date (LOD) to natural cooling than to warming, i.e. the delay of LOD caused by 1°C cooling is larger than the advance of LOD with 1°C warming. This might be because cooling leads to larger chilling accumulation and heating requirements for leaf onset, but this non-symmetric LOD response is partially offset by warming-related drying. Moreover, spring greening magnitude, in terms of satellite-based greenness and productivity, is more sensitive to LOD changes in the warming area than in the cooling. These results highlight the importance of considering non-symmetric responses of spring greening to warming and cooling when predicting vegetation-climate feedbacks.

Deeper topsoils enhance ecosystem productivity and climate resilience in arid regions, but not in humid regions.

Understanding the controlling mechanisms of soil properties on ecosystem productivity is essential for sustaining productivity and increasing resilience under a changing climate. Here we investigate the control of topsoil depth (e.g., A horizons) on long-term ecosystem productivity. We used nationwide observations (n = 2401) of topsoil depth and multiple scaled datasets of gross primary productivity (GPP) for five ecosystems (cropland, forest, grassland, pasture, shrubland) over 36 years (1986-2021) across the conterminous USA. The relationship between topsoil depth and GPP is primarily associated with water availability, which is particularly significant in arid regions under grassland, shrubland, and cropland (r = .37, .32, .15, respectively, p < .0001). For every 10 cm increase in topsoil depth, the GPP increased by 114 to 128 g C m-2 year-1 in arid regions (r = .33 and .45, p < .0001). Paired comparison of relatively shallow and deep topsoils while holding other variables (climate, vegetation, parent material, soil type) constant showed that the positive control of topsoil depth on GPP occurred primarily in cropland (0.73, confidence interval of 0.57-0.84) and shrubland (0.75, confidence interval of 0.40-0.94). The GPP difference between deep and shallow topsoils was small and not statistically significant. Despite the positive control of topsoil depth on productivity in arid regions, its contribution (coefficients: .09-.33) was similar to that of heat (coefficients: .06-.39) but less than that of water (coefficients: .07-.87). The resilience of ecosystem productivity to climate extremes varied in different ecosystems and climatic regions. Deeper topsoils increased stability and decreased the variability of GPP under climate extremes in most ecosystems, especially in shrubland and grassland. The conservation of topsoil in arid regions and improvements of soil depth representation and moisture-retention mechanisms are critical for carbon-sequestration ecosystem services under a changing climate. These findings and relationships should also be included in Earth system models.

Global water use efficiency saturation due to increased vapor pressure deficit.

The ratio of carbon assimilation to water evapotranspiration (ET) of an ecosystem, referred to as ecosystem water use efficiency (WUEeco), is widely expected to increase because of the rising atmospheric carbon dioxide concentration (Ca). However, little is known about the interactive effects of rising Ca and climate change on WUEeco. On the basis of upscaled estimates from machine learning methods and global FLUXNET observations, we show that global WUEeco has not risen since 2001 because of the asymmetric effects of an increased vapor pressure deficit (VPD), which depressed photosynthesis and enhanced ET. An undiminished ET trend indicates that rising temperature and VPD may play a more important role in regulating ET than declining stomatal conductance. Projected increases in VPD are predicted to affect the future coupling of the terrestrial carbon and water cycles.

New-generation geostationary satellite reveals widespread midday depression in dryland photosynthesis during 2020 western U.S. heatwave.

Emerging new-generation geostationary satellites have broadened the scope for studying the diurnal cycle of ecosystem functions. We exploit observations from the Geostationary Operational Environmental Satellite-R series to examine the effect of a severe U.S. heatwave in 2020 on the diurnal variations of ecosystem photosynthesis. We find divergent responses of photosynthesis to the heatwave across vegetation types and aridity gradients, with drylands exhibiting widespread midday and afternoon depression in photosynthesis. The diurnal centroid and peak time of dryland gross primary production (GPP) substantially shift toward earlier morning times, reflecting notable water and heat stress. Our geostationary satellite-based method outperforms traditional radiation-based upscaling methods from polar-orbiting satellite snapshots in estimating daily GPP and GPP loss during heatwaves. These findings underscore the potential of geostationary satellites for diurnal photosynthesis monitoring and highlight the necessity to consider the increased diurnal asymmetry in GPP under stress when evaluating carbon-climate interactions.

Temporal dynamics of ecosystem, inherent, and underlying water use efficiencies of forests, grasslands, and croplands and their responses to climate change.

BACKGROUND: Understanding temporal trends and varying responses of water use efficiency (WUE) to environmental changes of diverse ecosystems is key to predicting vegetation growth. WUE dynamics of major ecosystem types (e.g., forest, grassland and cropland) have been studied using various WUE definitions/metrics, but a comparative study on WUE dynamics and their driving forces among different ecosystem types using multiple WUE metrics is lacking. We used eddy covariance measurements for 42 FLUXNET2015 sites (396 site years) from 1997 to 2014, as well as three commonly used WUE metrics (i.e., ecosystem, inherent, and underlying WUE) to investigate the commonalities and differences in WUE trends and driving factors among deciduous broadleaf forests (DBFs), evergreen needleleaf forests (ENFs), grasslands, and croplands. RESULTS: Our results showed that the temporal trends of WUE were not statistically significant at 73.8% of the forest, grassland and cropland sites, and none of the three WUE metrics exhibited better performance than the others in quantifying WUE. Meanwhile, the trends observed for the three WUE metrics were not significantly different among forest, grassland and cropland ecosystems. In addition, WUE was mainly driven by atmospheric carbon dioxide concentration at sites with significant WUE trends, and by vapor pressure deficit (VPD) at sites without significant trends (except cropland). CONCLUSIONS: Our findings revealed the commonalities and differences in the application of three WUE metrics in disparate ecosystems, and further highlighted the important effect of VPD on WUE change.

Vegetation photosynthetic phenology dataset in northern terrestrial ecosystems.

Vegetation phenology can profoundly modulate the climate-biosphere interactions and thus plays a crucial role in regulating the terrestrial carbon cycle and the climate. However, most previous phenology studies rely on traditional vegetation indices, which are inadequate to characterize the seasonal activity of photosynthesis. Here, we generated an annual vegetation photosynthetic phenology dataset with a spatial resolution of 0.05 degrees from 2001 to 2020, using the latest gross primary productivity product based on solar-induced chlorophyll fluorescence (GOSIF-GPP). We combined smoothing splines with multiple change-point detection to retrieve the phenology metrics: start of the growing season (SOS), end of the growing season (EOS), and length of growing season (LOS) for terrestrial ecosystems above 30° N latitude (Northern Biomes). Our phenology product can be used to validate and develop phenology or carbon cycle models and monitor the climate change impacts on terrestrial ecosystems.

Large diurnal compensatory effects mitigate the response of Amazonian forests to atmospheric warming and drying.

Photosynthesis and evapotranspiration in Amazonian forests are major contributors to the global carbon and water cycles. However, their diurnal patterns and responses to atmospheric warming and drying at regional scale remain unclear, hindering the understanding of global carbon and water cycles. Here, we used proxies of photosynthesis and evapotranspiration from the International Space Station to reveal a strong depression of dry season afternoon photosynthesis (by 6.7 ± 2.4%) and evapotranspiration (by 6.1 ± 3.1%). Photosynthesis positively responds to vapor pressure deficit (VPD) in the morning, but negatively in the afternoon. Furthermore, we projected that the regionally depressed afternoon photosynthesis will be compensated by their increases in the morning in future dry seasons. These results shed new light on the complex interplay of climate with carbon and water fluxes in Amazonian forests and provide evidence on the emerging environmental constraints of primary productivity that may improve the robustness of future projections.

Persistent and enhanced carbon sequestration capacity of alpine grasslands on Earth's Third Pole.

The carbon sequestration capacity of alpine grasslands, composed of alpine meadows and steppes, in the Tibetan Plateau has an essential role in regulating the regional carbon cycle. However, inadequate understanding of its spatiotemporal dynamics and regulatory mechanisms restricts our ability to determine potential climate change impacts. We assessed the spatial and temporal patterns and mechanisms of the net ecosystem exchange (NEE) of carbon dioxide in the Tibetan Plateau. The carbon sequestration of the alpine grasslands ranged from 26.39 to 79.19 Tg C year-1 and had an increasing rate of 1.14 Tg C year-1 between 1982 and 2018. While alpine meadows were relatively strong carbon sinks, the semiarid and arid alpine steppes were nearly carbon neutral. Alpine meadow areas experienced strong increases in carbon sequestration mainly because of increasing temperatures, while alpine steppe areas had weak increases mainly due to increasing precipitation. Carbon sequestration capacity of alpine grasslands on the plateau has undergone persistent enhancement under a warmer and wetter climate.

Convergence and divergence emerging in climatic controls of polynomial trends for northern ecosystem productivity over 2000-2018.

Recent rapid warming has caused uneven impacts on the composition, structure, and functioning of northern ecosystems. It remains unknown how climatic drivers control linear and non-linear trends in ecosystem productivity. Based on a plant phenology index (PPI) product at a spatial resolution of 0.05° over 2000-2018, we used an automated polynomial fitting scheme to detect and characterize trend types (i.e., polynomial trends and no-trends) in the yearly-integrated PPI (PPIINT) for northern (> 30°N) ecosystems and their dependence on climatic drivers and ecosystem types. The averaged slope for the linear trends (p < 0.05) of PPIINT was positive across all the ecosystems, among which deciduous broadleaved forests and evergreen needle-leaved forests (ENF) showed the highest and lowest mean slopes, respectively. More than 50% of the pixels in ENF, arctic and boreal shrublands, and permanent wetlands (PW) had linear trends. A large fraction of PW also showed quadratic and cubic trends. These trend patterns agreed well with estimates of global vegetation productivity based on solar-induced chlorophyll fluorescence. Across all the biomes, PPIINT in pixels with linear trends showed lower mean values and higher partial correlation coefficients with temperature or precipitation than in pixels without linear trends. Overall, our study revealed the emergence of latitudinal convergence and divergence in climatic controls on the linear and non-linear trends of PPIINT, implying that northern shifts of vegetation and climate change may potentially increase the non-linear nature of climatic controls on ecosystem productivity. These results can improve our understanding and prediction of climate-induced changes in plant phenology and productivity and facilitate sustainable management of ecosystems by accounting for their resilience and vulnerability to future climate change.

Quantifying the impacts of land cover change on gross primary productivity globally.

Historically, humans have cleared many forests for agriculture. While this substantially reduced ecosystem carbon storage, the impacts of these land cover changes on terrestrial gross primary productivity (GPP) have not been adequately resolved yet. Here, we combine high-resolution datasets of satellite-derived GPP and environmental predictor variables to estimate the potential GPP of forests, grasslands, and croplands around the globe. With a mean GPP of 2.0 kg C m-2 yr-1 forests represent the most productive land cover on two thirds of the total area suitable for any of these land cover types, while grasslands and croplands on average reach 1.5 and 1.8 kg C m-2 yr-1, respectively. Combining our potential GPP maps with a historical land-use reconstruction indicates a 4.4% reduction in global GPP from agricultural expansion. This land-use-induced GPP reduction is amplified in some future scenarios as a result of ongoing deforestation (e.g., the large-scale bioenergy scenario SSP4-3.4) but partly reversed in other scenarios (e.g., the sustainability scenario SSP1-1.9) due to agricultural abandonment. Comparing our results to simulations from state-of-the-art Earth System Models, we find that all investigated models deviate substantially from our estimates and from each other. Our maps could be used as a benchmark to reduce this inconsistency, thereby improving projections of land-based climate mitigation potentials.

Optimal ranges of social-environmental drivers and their impacts on vegetation dynamics in Kazakhstan.

Kazakhstan is part of the Eurasian Steppes, the world's largest contiguous grassland system. Kazakh grassland systems are largely understudied despite being historically important for agropastoral practices. These grasslands are considered vulnerable to anthropogenic activities and climatic variability. Few studies have examined vegetation dynamics in Central Asia owing to the complex impacts of moisture, climatic and anthropogenic forcings. A comprehensive analysis of spatiotemporal changes of vegetation and its driving factors will help elucidate the causes of grassland degradation. Here, we investigated the individual and pairwise interactive influences of various social-environmental system (SES) drivers on greenness dynamics in Kazakhstan. We sought to examine whether there is a relationship between peak season greenness and its drivers - spring drought, preceding winter freeze-thaw cycles, percent snow cover and snow depth - for Kazakhstan during 2000-2016. As hypothesized, snow depth and spring drought accounted for 60 % and 52 % of the variance in the satellite-derived normalized difference vegetation index (NDVI) in Kazakhstan. The freeze-thaw process accounted for 50 % of NDVI variance across the country. In addition, continuous thawing during the winter increased vegetation greenness. We also found that moisture and topographic factors impacted NDVI more significantly than socioeconomic factors. However, the impacts of socioeconomic drivers on vegetation growth were amplified when they interacted with environmental drivers. Terrain slope and soil moisture had the highest q-values or power of determinant, accounting for ~70 % of the variance in NDVI across the country. Socioeconomic drivers, such as crop production (59 %), population density (48 %), and livestock density (26 %), had significant impacts on vegetation dynamics in Kazakhstan. We found that most of the pairwise interactive influences of the drivers exhibited bi-factor enhancement, and the interaction between soil moisture and elevation was the largest (q = 0.92). Our study revealed the optimal ranges and tipping points of SES drivers and quantified the impacts of various driving factors on NDVI. These findings can help us identify the factors causing grassland degradation and provide a scientific basis for ecological protection in semiarid regions.

Soil moisture regulates warming responses of autumn photosynthetic transition dates in subtropical forests.

Autumn phenology plays a key role in regulating the terrestrial carbon and water balance and their feedbacks to the climate. However, the mechanisms underlying autumn phenology are still poorly understood, especially in subtropical forests. In this study, we extracted the autumn photosynthetic transition dates (APTD) in subtropical China over the period 2003-2017 based on a global, fine-resolution solar-induced chlorophyll fluorescence (SIF) dataset (GOSIF) using four fitting methods, and then explored the temporal-spatial variations of APTD and its underlying mechanisms using partial correlation analysis and machine learning methods. We further predicted the APTD shifts under future climate warming conditions by applying process-based and machine learning-based models. We found that the APTD was significantly delayed, with an average rate of 7.7 days per decade, in subtropical China during 2003-2017. Both partial correlation analysis and machine learning methods revealed that soil moisture was the primary driver responsible for the APTD changes in southern subtropical monsoon evergreen forest (SEF) and middle subtropical evergreen forest (MEF), whereas solar radiation controlled the APTD variations in the northern evergreen-broadleaf deciduous mixed forest (NMF). Combining the effects of temperature, soil moisture and radiation, we found a significantly delayed trend in APTD during the 2030-2100 period, but the trend amplitude (0.8 days per decade) was much weaker than that over 2003-2017. In addition, we found that machine learning methods outperformed process-based models in projecting APTD. Our findings generate from different methods highlight that soil moisture is one of the key players in determining autumn photosynthetic phenological processes in subtropical forests. To comprehensively understand autumn phenological processes, in-situ manipulative experiments are urgently needed to quantify the contributions of different environmental and physiological factors in regulating plants' response to ongoing climate change.

The Grain-for-Green project offsets warming-induced soil organic carbon loss and increases soil carbon stock in Chinese Loess Plateau.

The dynamics of soil organic carbon (SOC) stock is a vital element affecting the climate, and ecological restoration is potentially an effective measure to mitigate climate change by enhancing vegetation and soil carbon stocks and thereby offsetting greenhouse gas emissions. The Grain-for-Green project (GFGP) implemented in Chinese Loess Plateau (LP) since 1999 is one of the largest ecological restoration projects in the world. However, the contributions of ecological restoration and climate change to ecosystem soil carbon sequestration are still unclear. In this study, we improved a soil carbon decomposition framework by optimizing the initial SOC stock based on full spatial simulation of SOC and incorporating the priming effect to investigate the SOC dynamics across the LP GFGP region from 1982 through 2017. Our results indicated that SOC stock in the GFGP region increased by 20.18 Tg C from 1982 through 2017. Most portion (15.83 Tg C) of the SOC increase was accumulated when the GFGP was initiated, with a SOC sink of 16.12 Tg C owing to revegetation restoration and a carbon loss of 0.29 Tg C due to warming during this period. The relationships between SOC and forest canopy height and investigations on the SOC dynamics after afforestation revealed that the accumulation rate of SOC could be as high as 24.68 g C m-2 yr-1 during the 70 years following afforestation, and that SOC could decline thereafter (-8.89 g C m-2 yr-1), which was mainly caused by warming. This study provides a new method for quantifying the contribution of ecological restoration to SOC changes, and also cautions the potential risk of LP SOC loss in the mature forest soil under future warming.

Decipher soil organic carbon dynamics and driving forces across China using machine learning.

The dynamics of soil organic carbon (SOC) play a critical role in modulating global warming. However, the long-term spatiotemporal changes of SOC at large scale, and the impacts of driving forces remain unclear. In this study, we investigated the dynamics of SOC in different soil layers across China through the1980s to 2010s using a machine learning approach and quantified the impacts of the key factors based on factorial simulation experiments.Our results showed that the latest (2000-2014) SOC stock in the first meter soil (SOC100 ) was 80.68 ± 3.49 Pg C, of which 42.6% was stored in the top 20 cm, sequestrating carbon with a rate of 30.80 ± 12.37 g C m-2  yr-1 since the 1980s. Our experiments focusing on the recent two periods (2000s and 2010s) revealed that climate change exerted the largest relative contributions to SOC dynamics in both layers and warming or drying can result in SOC loss. However, the influence of climate change weakened with soil depth, while the opposite for vegetation growth. Relationships between SOC and forest canopy height further confirmed this strengthened impact of vegetation with soil depth and highlighted the carbon sink function of deep soil in mature forest. Moreover, our estimates suggested that SOC dynamics in 71% of topsoil were controlled by climate change and its coupled influence with environmental variation (CE). Meanwhile, CE and the combined influence of climate change and vegetation growth dominated the SOC dynamics in 82.05% of the first meter soil. Additionally, the national cropland topsoil organic carbon increased with a rate of 23.6 ± 7.6 g C m-2  yr-1 since the 1980s, and the widely applied nitrogenous fertilizer was a key stimulus. Overall, our study extended the knowledge about the dynamics of SOC and deepened our understanding about the impacts of the primary factors.

Combining GOES-R and ECOSTRESS land surface temperature data to investigate diurnal variations of surface urban heat island.

The surface urban heat island (SUHI) phenomenon is characterized by both high spatial and temporal variability, while its diurnal (i.e., diel) variations have rarely been investigated because traditional satellites and sensors flying on polar orbits (e.g., Landsat, MODIS) have no diurnal sampling capability. Here we combined land surface temperature (LST) data from the Geostationary Operational Environmental Satellites (GOES-R) and the Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) to explore the diurnal variations of SUHI and thermal differentiation among various land covers over the Boston Metropolitan Area. With the combined use of the LST data from GOES-R and ECOSTRESS, we took advantage of the strengths of both GOES-R (i.e., high frequency in each day and night) and ECOSTRESS (i.e., much finer spatial resolution). The SUHI intensity of the urban-core and suburban areas both exhibited clear diurnal patterns for different seasons: a continuous increase in the SUHI intensity from sunrise to noon and a decrease thereafter to sunset, followed by a relatively low and constant intensity during nighttime. The LST contrasts among different land cover types were clearly larger in the daytime than at nighttime and peaked around midday. At noon in summer, the LST of 'Developed, High Intensity' was 2.6 °C higher than that of 'Developed, Medium Intensity', and about 4.6 °C higher than that of "Developed, Open Space" and "Developed, Low Intensity". Controlling the percent impervious surface in construction land at a relatively low level (e.g., below ~49%) could effectively alleviate the impacts of SUHI. Compared with GOES-R data, ECOSTRESS LST is suitable for monitoring the diurnal variations of intracity thermal environment at the subdistrict (or neighborhood) scale. Our study highlights the value of the combined use of geostationary satellite and ECOSTRESS LST in exploring the diurnal cycling of the SUHI, and can help inform urban planning and land-based climate mitigation policies in the context of climate change.